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Many shear flows follow a route to turbulence that has striking similarities to bifurcation scenarios 
in low-dimensional dynamical systems. Among the bifurcations that appear, crisis bifurcations are 
important because they cause global transitions between open and closed attractors, or indicate 
drastic increases in the range of the state space that is covered by the dynamics. We here study 
exterior and interior crisis bifurcations in direct numerical simulations of transitional plane Poiseuille 
flow in a mirror-symmetric subspace. We trace the state space dynamics from the appearance 
of the first three-dimensional exact coherent structures to the transition from an attractor to a 
chaotic saddle in an exterior crisis. For intermediate Reynolds numbers, the attractor undergoes 
several interior crises, in which new states appear and intermittent behavior can be observed. The 
bifurcations contribute to increasing the complexity of the dynamics and to a more dense coverage 
of state space. 


Numerical and experimental studies of pipe and plane 
Couette flow have demonstrated the significance of exact 
coherent structures and their bifurcations for the transi¬ 
tion to turbulence [1-4]. Typically, these states appear 
in saddle-node bifurcations and then undergo further bi¬ 
furcations. Initially, most of their complexity lies in the 
temporal dynamics, so that they are better character¬ 
ized as chaotic rather than turbulent. With increasing 
Reynolds number, more temporal and spatial degrees of 
freedom are activated, until the complexity of a turbu¬ 
lent flow is established. Parallel to the increase in com¬ 
plexity comes a growth of the parts of state space that 
participate in the chaotic and turbulent dynamics. Stud¬ 
ies of low-dimensional dynamical systems have revealed 
many routes to this increased complexity [5-8]. Several 
of them have already been discussed in the context of 
high-dimensional fluid systems, e.g. in the cases of plane 
Couette flow [2] or pipe flow [4, 9]. One contribution of 
the present study is to document similar phenomenology 
in another canonical fluid system, plane Poiseuille flow 
(PPF). A second one is the demonstration of interior cri¬ 
sis and their contribution to increasing the complexity of 
the attractor and of the state space region covered by it. 

PPF is the pressure driven flow between two parallel 
plates and differs from plane Couette flow and pipe flow 
because of the presence of a linear instability to trans¬ 
verse vortices, the so-called Tollmien-Schlichting modes 
[10-12]. It occurs at a critical Reynolds number of 
5772.22 for a streamwise wavenumber a of 1.02056 (based 
on the center-line velocity and half the gap width), as de¬ 
termined by Orszag [13]. The bifurcation is subcritical, 
and reaches down to about Re ~ 2700 [14, 15] (for differ¬ 
ent wavelength). However, several experiments and nu¬ 
merical simulations show that turbulence occurs already 
at Reynolds numbers around 1000 [16-18], and hence well 
below the onset of Tollmien-Schlichting modes. Thus, the 
linear instability cannot explain the observed turbulence 
at low Reynolds numbers and the situation becomes anal¬ 
ogous to that in plane Couette and pipe flow. 

In order to determine the relevant saddle-node bifur¬ 
cation in PPF we use the method of edge tracking, as 


described in [19], see also [20]. The method traces the 
time-evolution of initial conditions and uses bisection be¬ 
tween an initial condition that returns to the laminar 
profile and one that becomes turbulent to approximate 
one on the laminar-turbulent interface. In most cases 
the state evolves towards a simple attractor, such as a 
fixed point or a simple periodic orbit. It is then possible 
to continue the edge state in Reynolds number around 
the saddle-node bifurcation and to identify the upper 
branch solution. In recent work for plane Couette [2] 
and pipe flow [4] it was shown that the upper branch 
of the edge state undergoes various bifurcations result¬ 
ing in a chaotic attractor. A boundary or exterior cri¬ 
sis ultimately destroys this stable attractor and creates 
the observed transient turbulence with its characteristic 
exponentially distributed lifetimes. The observed phe¬ 
nomenology is similar to what has been described and 
discussed in the context of chaotic dynamical systems 
[7, 8]. 

We will here show that this scenario is also present in 
plane Poiseuille flow and that another type of crisis bifur¬ 
cation, the interior crisis, provides a mechanism by which 
the part of state space occupied by the chaotic attractor 
can increase. Furthermore, we will discuss mechanisms 
for the observed increase of lifetimes of chaotic transients 
[ 21 ]- 

For our numerical simulations we use the Channelfow- 
code [22], The Reynolds number Re = Uod/u for the 
system is based on half the distance between the plates 
d , the maximum velocity of the laminar profile Uq and 
the kinematic viscosity v. We take a coordinate system in 
which x points in the streamwise, z in the spanwise and 
y in the wall-normal direction. With the above choices 
for the dimensionless units the laminar profile becomes 
ui(y) = (U(y), 0,0) with U(y) = 1 — y 2 . The total flow 
field ut can be written as the sum of the laminar profile 
and a fluctuating component, Ut = ui + u. All simula¬ 
tions in the paper are performed for constant mass flux 
and with no-slip boundary conditions at the walls. The 
calculations are restricted to a computational domain of 
length 2-7T, width n (and height 2) in a subspace that is 
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FIG. 1. Bifurcation diagram for the travelling wave TWe V z- 
Solid lines indicate a stable state, dashed lines an unstable 
state. Chaotic and periodic states are indicated by clouds 
of points obtained by plotting minima and maxima of their 
energy densities in the course of time. The red dots in¬ 
dicate the bifurcation point of TWe V z and PO yz , respec¬ 
tively. The insets show the average streamwise velocity in 
the spanwise-wallnormal plane for the lower and the upper 
branch of TWe V z at Re = 830. The colors indicate low (blue) 
and high (red) velocity regions. 

symmetric to reflections at the midplane and to spanwise 
reflections at the plane defined by z = 0: 

s v : [u, v, w\{x , y, z) = [u, -v, w]{x, -y, z) (1) 

s z : [u, v, w] (x, y, z) = [u, v, - w] (x, y, -z) (2) 

As in other studies [2, 4] the restriction to a symmetric 
subspace stabilizes the exact coherent structures. 

The numerical resolution is N x x N y x N z = 48 x 65 x 48 
modes. We checked higher values of N x and N z and 
found no significant changes in the bifurcations. The 
exact coherent structures and the location of the bifur¬ 
cations vary with domain size, but the phenomenology 
is similar. The same applies to the localized coherent 
states and their bifurcations that appear in wider do¬ 
mains. Especially, in large domains the states become 
localized [23, 24]. 

Using the technique of edge tracking it is possible to 
identify the edge state [19, 25] of this system. A trajec¬ 
tory on the laminar-turbulent boundary quickly reaches 
a states of constant energy. Since a stationary state is 
ruled out on account of the non-zero mean flow, the at¬ 
tractor in the laminar-turbulent boundary is a travelling 
wave. Indeed, a Newton search [26] for a relative equilib¬ 
rium converges to a traveling wave, henceforth referred 
to as TWe V z■ The form of the state is indicated in the 
inset in Fig. 1. The travelling wave has the same symme¬ 
tries as the mirror-symmetric travelling wave previously 
described by Nagata and Deguchi [271 and Gibson and 
Brandt [28]. 

We use a continuation method (see e.g. [29]) to follow 
the solutions in Reynolds number around the saddle-node 


FIG. 2. Magnification of the attractor in Fig. 1 in order to 
highlight the interior crisis (blue line) and the exterior crisis 
(red line). The states are visualized by by plotting minima of 
E(u) along trajectories. 

bifurcation at R esN = 641. A stability analysis of the 
travelling wave shows that in the symmetry subspace the 
lower branch has one unstable eigenvalue and the upper 
branch is stable for 641 < Re < 707. 

At Re = 707 the upper branch undergoes a Hopf bifur¬ 
cation that creates a stable relative periodic orbit PO yz . 
This orbit undergoes a Neimark-Sacker bifurcation at 
Re = 761.5 that creates a stable torus. In further bi¬ 
furcations a chaotic attractor is generated. By plotting 
minima and maxima of the energy 

i f Lz r 1 r Lx 

E(u) = / / / u 2 dxdydz (3) 

ZL/x-L/z Jo J — 1 Jo 

of a trajectory on the attractor we are able to map out the 
bifurcation diagrams also in chaotic regions, as shown in 
Fig. 1. The mapping of the chaotic region becomes fea¬ 
sible due to the restriction by the shift-and-reflect sym¬ 
metry that stabilizes the states. In the full system the 
entire bifurcation structure persists within an unstable 
subspace. 

The magnification of the chaotic attractor in Fig. 2 
highlights the two phenomena we want to focus on here: 
Slightly above Re = 785 (blue line) the size of the at¬ 
tractor expands and covers a larger fraction of the inter¬ 
val, and slightly below Re = 786 (red line) it disappears. 
Both changes are connected with crisis bifurcations [7, 8] , 
an interior crisis in the first case, and an exterior crisis 
in the second case. 

Slightly above Re = 785 the points on the attractor 
suddenly spread over a wider region, covering the area in 
state space with an energy E between 0.023 and 0.044. 
However, these parts of the state space are only visited 
occasionally, so that the points are less dense than in 
other parts. The reason for the sudden enlargement is a 
so-called interior crisis bifurcation [30, 31], where a new 
state appears and new links to the attractor form. The 
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appearance of the new states can be seen in the time 
series in Fig. 3. Just before the crisis, the range of the 
trajectories is limited to the interval [0.033,0.042]. 
Slightly above the crisis, excursions to lower values occur, 
with their number increasing with Re . 

The type of transition can be determined from the dis¬ 
tribution of times spent in the different regions. The 
state space region covered by the attractor before the 
crisis, referred to as phase A, contains trajectories that 
never drop below a threshold in energy, here taken to 
be E t = 0.031. If the trajectory drops below E t , tra¬ 
jectories enter a different phase B, occupying a different 
region in state space. An indicator for phase B are the 
repeated excursions to values below E t . Accordingly, if 
no excursions are noted for more 250 time units, we con¬ 
clude that the system has returned to phase A. With 
this prescription one can determine the distribution of 
times in phase A as shown in Fig. 4a. The plot contains 
data from trajectories with a total length of 5 • 10 6 time 
units. The data are shown semi-logarithmically, so that 
the times are compatible with an exponential distribu¬ 
tion, as expected for an interior crisis [32, 33]. We then 
fit an exponential decay to the distribution to obtain the 
characteristic trapping times t a in phase A and plot them 
versus Reynolds number in Fig. 4b. Approaching the cri¬ 
sis point from above, the time in phase A diverges since 
B is never visited. According to [ 8 , 32] the characteristic 
time varies as 

t a oc (Re - Re/c) -7 (4) 

with an exponent 7 . We use Re/c = 785.1, as it is the 
lowest values of Re for which we observe excursions to 
phase B , and fit the the exponent to the data. In the 
present case we obtain a good fit to the data with 7 = 0 . 8 . 

The exponents for the interior crisis (as well as those 
of the exterior crisis) are expected to be larger or equal 
to 1/2 for a smooth dynamical system [ 8 , 34]. For one 
dimensional maps with a quadratic maximum they can 
be shown to be exactly 1/2 but in higher dimensional 
systems the folding of the manifolds contributes to the 
dimensions, and higher exponents have been found [7, 8 , 
35, 36]. Since the exponents depend on the eigenvalues at 
the point of bifurcation different crises can show different 
exponents, even in the same physical system. 

Typical turbulent trajectories show an enormous tem¬ 
poral and spatial complexity that is difficult to create in a 
sequence of simple Neimark-Sacker or period-doubling bi¬ 
furcations. As is evident form Fig. 3 the dynamics of the 
system is rather regular (but not periodic!) before the in¬ 
terior crises and becomes increasingly more complicated 
(both in the range covered and in the complexity of the 
time-signal) as the Reynolds number increases. Thus, the 
interior crisis bifurcation increases the complexity of the 
chaotic trajectories more dramatically than other local 
bifurcations and are an important contribution towards 
more turbulent time evolutions. 

The second phenomenon we want to address here is 
the change in the dynamics near Re = 786, where the 
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FIG. 3. Energy traces E(t ) for trajectories near the interior 
crisis, (a) Re = 785.0, slightly below the crisis, (b), (c), and 
(d) are for Re = 785.46, 785.75, and 785.9, respectively, above 
the crisis. They show the characteristic intermittent bursts. 

chaotic attractor suddenly disappears. Here, the attrac¬ 
tor collides with the lower-branch state and turns into a 
chaotic saddle in a boundary or exterior crisis bifurcation 
[30] . It is a generic property of a chaotic saddle that the 
survival probabilities are exponentially distributed. To 
quantify this defining property of the boundary crisis, 
the survival probabilities for Re > Rexc are calculated 
using the methods described by [37]. The survival prob¬ 
abilities are clearly exponential distributed with charac¬ 
teristic lifetimes that depend on the Reynolds number, 
as shown in Fig. 5. As in the case of the interior crisis 
they diverge as 

r oc (Re — Rexc ) _<5 (5) 

for Re near the Reynolds number Rexc of the crisis bi¬ 
furcation. We fix Rexc = 785.95, since this is the lowest 
Reynolds number where we observed a trajectory that de¬ 
cays after showing transient chaotic dynamics for a long 
time. Best fits to the data are obtained for <5 = 1.5, 
as expected for an exterior crisis. A dense sampling of 
initial conditions in the state space of the system com¬ 
bined with a fine scan of Reynolds numbers in the range 
between Re = 778.3 and Re = 780.6 reveals a small at¬ 
tractor Ai a inside of A±. This attractor disappears at 
Re = 780.6 in another boundary crisis bifurcation and 
above this Reynolds number initial conditions exist that 
transiently visit the saddle created by the boundary crisis 
of the attractor An, before suddenly settling down on A±. 
A example for such a trajectory is shown in Fig. 6 . Fur¬ 
thermore, no initial conditions can be found that tran¬ 
siently visit the saddle before becoming laminar. This 
behavior is strong evidence that A la lies completely in¬ 
side of the basin of A/. 

The presence of the chaotic saddle created in the 
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FIG. 4. Intermittency near the interior crisis, (a) Probabil¬ 
ity to stay t time units in the pre-crisis phase A. The times 
are exponentially distributed and the characteristic timescale 
ta increases with decreasing distance to Re/c- (b) Variation 
of characteristic times with Reynolds number. The continu¬ 
ous line shows an algebraic fit to equation (4) with 7 = 0.8. 
The inset shows the data on a doubly logarithmic scale where 
on the abscissa the distance to the Reynolds number of the 
interior crisis is used. 



Re 


FIG. 5. Characteristic lifetimes r vs. Reynolds number above 
the boundary crisis. The continuous line is an algebraic fit to 
equation (4) with <5 = 1.5. The inset shows r for a larger 
range in Re. The points are connected to guide the eye only. 
Regions where a stable attractor exists are shaded grey; the 
lifetime r is infinite in these regions. 



FIG. 6 . Transition between saddle A\b and attractor Ai in 
the energy trace of a trajectory at Re = 780.7. The trajectory 
stays on the chaotic saddle du for about 18000 time units 
before it suddenly switches to the stable attractor Ai. 


boundary crisis of A\ a should lead to a second slope in 
the lifetime distribution as also seen in [36]. But since 
the basin of A\ a is very small compared to A\ this slope 
does not influence the characteristic lifetimes in Fig. 5. 
The lifetimes for a larger range in Re are shown in the in¬ 
set of Fig. 5. Following the boundary crisis the lifetimes 
of the chaotic transients first decrease, then start to in¬ 
crease again around Re « 815, and diverge at Re = 828, 
where a second stable attractor (A 2 ) appears. At slightly 
higher Re another attractor (A 3 ) appears so that includ¬ 
ing the laminar state for a small range in Re the system 
has three attracting states. A 3 and A 3 disappear in a 
boundary crises at Re = 837.5 and Re = 841.8, respec¬ 
tively. After the boundary crisis of A 3 the lifetimes drop 
to even lower values than before the appearance of A 2 . 
They decrease until Re = 930, where a lifetime of 126 is 
reached. Afterwards lifetimes increase again and eventu¬ 
ally diverge at Re = 1087, where another attractor A 4 
appears. The attractors A 2 -A 4 appear in regions of the 
state space occupied by the large saddle created in the 
boundary crisis of A 3 , as was checked using slices of the 
state space as in [36]. 

The crisis bifurcations analyzed here for PPF extend 
previous observations on Couette flow [36] and also pipe 
flow [38] in that they provide further examples of smaller 
chaotic saddles inside larger outer saddles and of bifur¬ 
cations that contribute to an increase of the character¬ 
istic lifetimes and eventually to a more complex tempo¬ 
ral dynamics. Moreover, the interior crises contribute to 
a more dense coverage of the state space of the system 
by the dynamics. Thereby, they pave the way for the 
transition to a chaotic saddle when the attractor collides 
with the saddle from the original saddle-node bifurcation. 
The connection to the phenomenology of low-dimensional 
dynamical systems and the appearance in a number of 
canonical flows suggests that this transition scenario is 
typical for the transition in shear flows. 
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